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A general self-consistency approach allows a thorough treatment of the corrections to the mean- 
field approximation (MFA). The natural extension of standard MFA with the help of a cumulant 
expansion leads to a new point of view on the effective field theories. The proposed approach can 
be used for a systematic treatment of fluctuation effects of various length scales and, perhaps, for 
the development of a new coarse graining procedure. We outline and justify our method by some 
preliminary calculations. Results are given for the critical temperature and the Landau parameters 
of the </!> 4 -theory - the field counterpart of the Ising model. An important unresolved problem of the 
modern theory of phase transitions - the problem for the calculation of the true critical temperature, 
is considered within the framework of the present approach. A comprehensive description of the 
ground state properties of many-body systems is also demonstrated. 
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I. INTRODUCTION 

This investigation is focused on the correspondence be- 
tween microscopic models of phase transitions and their 
quasi-macroscopic (field-theoretic) counterparts. Here 
we shall outline a new self-consistency approach to a more 
accurate derivation of effective field theories from micro- 
scopic models defined on lattices. 

Our method is general and can be used for a wide class 
of microscopic models but for a concreteness here we shall 
illustrate our approach with the Ising model (IM) , given 
by 

1 N 

^( s ) = -g J v s i s i> (!) 

ij 

where s = {s,} denotes a lattice "field", Si = ±1, and 
the interaction constant of ferromagnetic type Jy = 
J{V ~ il) > depends on the intersite distance \i — j\ 
in a regular _D-dimensional lattice of N sites ("spins" or 
pseudo-spins). Note, that J(0) = Ju = 0. 

We shall follow the main path of the phase transi- 
tion theory, where the effective (quasi-macroscopic) field 
Hamiltonians (alias Ginzburg-Landau (GL) free energies) 
are derived with the help of two systematic methods: (i) 
Hubbard- Stratonovich transformations (HST) and, (ii) a 
mean- field (MF) like procedure [3, 0, Q . Here we pro- 
pose a more thorough approach based on a convenient 
generalization of (ii). 

The known field theories exhibit both success and fail- 
ure in the description of many-body systems defined by 
microscopic models. For example, we believe that the 
renormalisation group methods of the modern theory 
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of phase transitions |l| based on the </> -theory yield a 
quite convenient description of the scaling and univer- 
sality properties of IM but we cannot be certain that 
the field theory satisfactory describes important non- 
universal properties, such as, for example, the critical 
temperature T c and the lower critical dimensionality Dr, 
(for IM, Dr, = 1, whereas within the (/> 4 -theory, Dl = 2). 
In fact the field theory fails along this line of studies. 
For example, the fluctuation shift (AT C ) / of T c predicted 
within the one-loop approximation for the </> 4 -theory is 
a very small and, hence, unrealistic, while in the higher 
orders of the loop expansion this shift turns out infinite 
and its calculation needs a special renormalization. As a 
result, the problem for the value of T c within the frame- 
work of the present field theory of phase transitions re- 
mains unresolved. Here we shall show the genesis of this 
problem and present a satisfactory solution. 

The reason for the mentioned difficulties of the field 
theory approach is in the quite simplified coarse-graining 
procedures used in the derivation of the GL effective 
Hamiltonians from microscopic models. HST is ap- 
plied together with the long-wavelength approximation 
(LWLA), namely, (kao) <C 7r, where k = \k\ is the mag- 
nitude of the wave vector k = (ki, &d), ao is either 
the lattice constant or, generally, the mean inter-particlc 
distance. LWLA leads to a correct expression for the 
Ornstein-Zernicke correlation function but cuts the short- 
range (high-energy: e(k) ~ fc 2 ,A < k < (j)/do) interpar- 
ticle correlations of fluctuation type that have the main 
contribution to the shift (AT C )/; A -C Tr/ao is the upper 
cutoff for k within LWLA. 

Within the approach (ii) we usually say that we ne- 
glect the fluctuations of the physical quantities from their 
equilibrium values. This is not entirely true. It seems im- 
portant to emphasize that the latter are values in MF ap- 
proximation (MFA) and, hence, they are incorrect. Thus 
the fluctuations are defined towards incorrect statistical 
averages, and their contribution to the free energy of the 
system cannot be accepted as an entire fluctuation effect. 
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The latter can be correctly evaluated, if we are able to 
define the fluctuations as variations towards exact statis- 
tical averages calculated by the Hamiltonian (TTJ). This 
task seems unsolvable, but the present paper makes a 
step of improvement of the theory along the same direc- 
tion. 

Note, that the methods (i) and (ii), no matter of the 
difference between them, lead to the same GL effective 
field theory. Both methods use LWLA. For the method 
(ii) LWLA seems to me obligatory because of the fol- 
lowing important argument. HST can be applied only to 
positively definite matrices (Jij) but this condition is not 
satisfied by quite important lattice models, such as the 
nearest-neighbour (raj) 1M. In LWLA, however, the same 
interaction matrix Jij ;nra of IM is modified to a form that 
is positively definite. Besides, in both methods, LWLA 
is used to help the derivation of the quasi-macroscopic 
(fluctuation) and macroscopic (thermodynamic) proper- 
ties, i.e., the LWLA is a tool of a "coarse-graining" pro- 
cedure for the many-body system. But in both cases (i) 
and (ii), the micro- and mesoscopic interparticle correla- 
tions are ignored. Here we shall use a generalization of 
the method (ii) in order to improve this disadvantage of 
the theory. In this way we shall present more accurate 
calculation of the LG parameters (vertices) of the effec- 
tive field theory of Ising systems. Besides we shall show 
explicitly for the first time the mechanism of "statistical 
correlation" which leads from the two-site (i — j) triv- 
ial correlation presented by the initial interaction Jy- to 
short-, meso-, and large-scale effective multiparticle cor- 
relations of fluctuation type - fluctuation correlations. 
Thus we shall establish and develop for the first time a 
new "coarse-graining" procedure, and this is the main 
aim of our report. Besides, we shall show that the GL 
parameters of the 4 -theory acquire (l/z)-corrections for 
both short-range and long-range interactions J(\i— j\) [H- 
These corrections will be presented to second order in 
(1/zf. 

The I /z-expansion has been introduced by R. Brout [J] 
and applied in calculations of T c and thermodynamic sus- 
ceptibilities [lEHBSIISdEd- We are not aware 
of another relevant work along this line of research ex- 
cept for our recent investigation 0, H[ ; the latter will be 
used in our investigation. Remember, that the (Brout) 
approach is a development of an older method - the Kirk- 
wood method of cumulants (semi- invariants) [lij : hence, 
we shall follow a cumulant expansion which is well known. 

We shall generalize the Brout approach in a way that 
makes possible to derive effective field theories, and this 
will allow us to reveal new and surprising features of 
many-body systems. In the Brout scheme, the mean 
("molecular") field is spatially (i-) independent. We find 
that there are no physical reasons for this assumption and 
extend the "mean-field concept." Within our approach, 
the so-called mean field is spatially dependent up to the 
moment when one should find the actual ground state. 
Then the fluctuation phenomena occur with respect to 
this ground state. In our approach the GL parameters 



and, hence, the ground state have a more precise evalu- 
ation. 

Our consideration has been performed for interactions 
J(\i — j\) of a quite general type, namely for all interac- 
tions that can be presented by the equality 

Jij = zJo = J, (2) 

3 

where Jo is an effective exchange constant, and z is an 
effective "coordination" number (number of interacting 
neighbors). In particular cases our results will be referred 
to the most common case of nn interactions (then n = 2D 
for simple lattices). For a simplicity, here we shall as- 
sume that IM is defined by a simple cubic (sc) lattice. 
The consideration can be easily expended to other types 
of regular lattices, as well as to irregular lattices with 
certain forms of quenched disorder, for example, random 
potential [l[. We introduce the interaction radius by 
Rint ~ z x l B - a quantity, which is equal to the so-called 
zero-temperature correlation length (see Section 2). Let 
emphasize that our results can be rederived without diffi- 
culties in the presence of an external field {hi} conjugate 
to the lattice field {si}. 

Other theories intended to improve the coarse graining 
procedures of many-body systems are, for example, the 
hierarchical reference theory [l5| and a method of collec- 
tive variables The former is based on a particular 
renormalization procedure and yields interesting results 
for the effective critical exponents of continuous phase 
transitions. The latter theory uses a variant of cumulant 
expansion and collective variables which, in certain limit, 
form the order parameter field. In these features the col- 
lective variable theory [l6| resembles our approach but, 
generally, in contrast to the latter, it significantly departs 
from the original Brout idea. The mentioned theories 
does not use 1 / z-expansion and for this reason and apart 
from some common aims and general ideas, cannot eas- 
ily be compared with the present approach, which is a 
very direct extension of the many-body theory beyond 
the standard mean- field ( "tree" ) approximation. 

In Sec. II we present a general approach to the treat- 
ment of fluctuation correlations at various length scales 
and discuss aspects of the usual theory that corresponds 
to the lowest order approximation of a perturbation ex- 
pansion of cumulant type. In Sec. Ill we investigate the 
higher orders of the mentioned perturbation expansion 
and demonstrate several new features of the effective field 
theory. Our main results are summarized and discussed 
in Sec. III.D - III.G. 



II. PRESENT STATUS OF THE EFFECTIVE 
FIELD THEORY 

A. General scheme 

The equilibrium free energy of IM as a function of the 
temperature T and the field configuration h = {hi} is 
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given by 



G(T) = -/r 1 ln{Tre-' w ^}, 



(3) 



where f3~ 1 = ksT, and the Trace is over the allowed 
lattice configurations s = {si}. Note, that the equilib- 
rium values of the physical quantities are calculated as 
averages with respect to the statistical ensemble based 
on the Hamiltonian ([T]) and, in particular, the statis- 
tical averages of type (si...Sj) are obtained as deriva- 
tives of the partition sum whereas the irreducible aver- 
ages ((si...Sj)) = ((sj — (si))...(sj — (sj))) are obtained as 
derivatives of the Gibbs free energy ((3|). 
Let us introduce the shift 

Si =<f>i + Ssi, (4) 

where the lattice field (pi is an arbitrary (auxiliary) field 
configuration that is not necessarily associated with the 
averaged spin (sj) at site i, and the "fluctuation" 8si 
is merely the difference (sj — (pi). The identification of 
(pi with a statistical average (...) over the full Hamilto- 
nian ([1]), or, with a statistical average (...)o corresponding 
to another ensemble as well as the interpretation of 5si 
as a fluctuation around (si) (or (sj)o) may be a matter 
of further considerations. At this stage 4>i and Ssi are 
auxiliary variables which obey ^ and are not referred 
to concrete physical quantities. 

Following a standard procedure (see, e.g., Refs. [l], 0) 
we obtain the following effective non-equilibrium free en- 
ergy 



H{(P)=H {(P)+Hf{(P), 

with 4> = {<£»}, 

W) = \ E hiMi - ^ E ln t 2ch ^ c 

ij 



Here 



E^ 



(5) 



(6) 



(7) 



is the "mean" (molecular) field, and 



Hf(4>) = -/TMn^exp 



f3 x 



is the "fluctuation" part. As usual, we shall often call 
the free energy an "effective Hamiltonian" (see, also, 
Ref. [13). 

In ([5]), (...)o denotes a statistical average over an en- 
semble defined by the auxiliary ("MF") Hamiltonian 



H a {(p,s) = -E° lSi ' 



(9) 



the respective partition function and (nonequilibrium) 
free energy are given by Z a {4>) = Tr [exp(— PTi. a )], and 



have 



-j3 1 lnZ a ((p). For this simple ensemble, we 



»<}o =th [I3a % {ct>)]. 



(10) 



The calculation of averages of type (si...sj)o as well as "ir- 
reducible" averages of type ((Ssi...Ssj))o is also straight- 
forward. These averages represent a form of fluctuation 
correlations, but they are just an auxiliary theoretical 
tool rather than real objects. In contrast, the real ob- 
jects, namely, (full) statistical averages (...) and ((...)) 
within the total Hamiltonian ([T]) cannot be exactly cal- 
culated. 



B. Usual theory 

In the framework of the usual theory, the "fluctuation" 
term Hf is ignored. This is the MFA. In the present 
general format of the theory, we have N (self-consistency) 
equations of state (&H/d(pi) — at fixed T (and h = 
{hi}) — 0) - one equation per a lattice vertex i: 



£jy{^-th[/3 Oi $)]} = 0, 



(11) 



where <f> = {(pi}. It is easy to see that the number (z — 1) 
of nonzero terms (i ^ j) in all ./V sums is equal to 
the number of nonzero interaction constants acting on the 
site i: J(\i-j\) > for a < < Rinu J(\i~j\ = 

for \i — j\ = R > Rint- The "equations of state" (TTTj) can 
be written in the simple form 



= th \j3cti 



(12) 



The equivalence of pT|) and (TT2]) can be easily proven for 
any number N > 1. 

From (fini) and (TT2]) we obtain 



) = th [{3ai{4>)] 



(13) 



Thus in this quite general form of MFA (Hf ~ 0), we 
have: { Sl ) = fa (...) = (...) , and G(T,h) = H $) 
is the equilibrium free energy that corresponds to ex- 
trema (including minima) <p of H((p) ~ Tio{(p) ~ the non- 
equilibrium MF free energy given by the non-equilibrium 
(arbitrary) order parameter field (pi. By (pi from (jlip we 
denote the equilibrium configuration of the latter; here- 
after the "bar" of the equilibrium value (pi of (pi will be 
often omitted. Now one may perform a Landau expan- 
sion for small (pi in order to get other known forms of the 
MF theory. 

However, at the present stage of consideration we are 
interested in some more generality and for this reason 
we continue our discussion of the non-equilibrium free 
energy functional 



G(T/4) = H(T/(P) » H Q {T/(p) 



(14) 
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as given by Eq. ([6]). In (fl"4"f . the dependence of the non- 
equilibrium free energy G on <fi is denoted by "/</>" be- 
cause of the more special role of this variable, namely, 
the variation 8G should be zero at thermal equilibrium 
and from this condition one obtains the possible thermal 



equilibria <f> (alias "self-consistency condition"). 

The expansion of the log— term in ([6|) up to fourth or- 
der in <j)i yields the known result for the lattice <pf— theory 
of the Ising model: 



H (<P) 



- ^ Jij Jik<f>j< 
ijk 



r 



12 ^ lJ 

ijklm 



JiiJikJilJim ( l ) j ( t ) k<j>l< l 



(15) 



One may apply LWLA to this form of the theory but we 
shall follow a different path. 



C. Continuum limit 

Using the rule 

J2h = p[ d ° X f( X ) = pj dx f( X ) > ( 16 ) 

where p = (N/V) we can write Eq. © in the form 

J = p J d D RJ(R), (17) 

with R — \R\; R — (x — y). Note, that the factor p in 
(fT9|) - (|21|) can be avoided and this is the usual practice. 
In the latter case the physical dimension of the respective 
physical quantity is changed by a factor [V] ~ [L] D \ for 
example, = [/(a;)]. Of course, one may use both 

variants. 

In the continuum limit the effective Hamiltonian ^ 
corresponding to a zero external field (h — 0) takes the 
form 



Uo = I d u x6(x)I\(h(x) 



(18) 



where 



p/T 1 / d D x\nch{(3- 1 I[(j)(x)}} 



I[<j>(x)]=p J d D yJ(R)4(y) 



(19) 



Now we apply LWLA in the form 

[#r)-0(y)] 2 « \<t>(x)ct>(y)\ (20) 
and under this assumption truncate the Taylor expansion 

'(y) = d?{x) + Y d 4^R a (21) 



dx a 

+ 1 y d 2 cj>{x) 
^ x dXadx/3 



a=l 

D 



R a Rt 



to the second order in R = {R a }. 

Using the approximation (|20|) . (fl9|) becomes 



J 

2D 



V 2 c 



I[<P(x)} = J4>(x) 
where J is given by pT|) . and 

J = p { d D RJ(R)R 2 



(22) 



(23) 



Now we substitute (|2"2"j) in (18), perform the expansion 
up to order <t) 4 (x) and to second order in V(j)(x). Besides, 
we should keep in mind, that in expansion in powers of 
<pi we cannot distinguish between T and T c0 except for 
the (f)f— term where the difference between T and T c q 
should be kept only to the lowest nonvanishing order; in 
our case, this is the first order in (T — T c q): see, e.g., 
Ref. [l|. Following these notes, we perform at a certain 
stage of the calculation an integration by parts with the 
convenient boundary condition \J(j){x) — and obtain 
the well known GL effective Hamiltonian 



Ho = p 



d D , 



|[V0(a;)] 2 + ^cj> 2 (x) +u Q ^(x) 



with 



Co 



— J, ro(T) 
2D ' y ' 



k B (T-T cl 



u 



J_ 

12' 



(24) 



(25) 



Here T c q — (J/ks) and terms of order to = (T — 
T c o)/T c q -C 1 have been neglected in Co and ur> 0, i.e. 
these two parameters are calculated at T c o [3, Q • 

To clarify the result (HH)) - (|23|) we shall mention that 
J(R) for R > Ri n t is very small and can be ignored. 
Setting J{R) ~ Jo in dHj) ; comparing the result with 
J = zJ from (2J), and noticing that p = 1/v 



l o ■> 



one 



obtains z ~ (Ri n t/ao) D as should be, and J w JoR? 



In the same way one gets J w JR 2 nt ■ 



The results (|24|) - (f25|) show that the energy of the 
spatially dependent configurations of the field depends 
on the interaction radius. The latter serves as a coher- 
ence (correlation) length of the field (j>(x); see the pa- 
rameter Co given in (|25p . In order to clarify this point, 
let us consider the so-called zero-temperature correlation 
length 0, defined by £ = £(T = 0) - hc /r (0)] 1/2 . 
Using (|2"5"|) and T c0 = J/ks we obtain ^ = RintV^D- 
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D. Discussion 

Note that the temperature range of validity of these 
considerations is to(T) <§; 1 and (ka ) <C w. The spa- 
tially dependent fluctuations correspond to a higher en- 
ergy than the uniform configuration, and hence the latter 
contains the deepest (global) minima cf> of the effective 
free energy nevertheless we have written N "equations 
of state" as a result of the minimization of the effective 
Hamiltonian. Now one can easily show that the varia- 
tion of the Hamiltonian (12411 with respect to the field will 
give again spatially dependent solutions but in the usual 
theory they are interpreted as "spatially dependent fluc- 
tuations" which, together with the uniform fluctuation 
5(f> = ((f) — <j>) towards the stable state <j>, are all fluctu- 
ations in the system in LWLA. But we know that this 
picture contains the approximation (sj) ~ (sj)o = 4>- 

Our point of view is the following. The averages ()o 
should not be taken very seriously. They are an auxil- 
iary tool in our consideration, and are not the final aim of 
our investigation. We have used these averages only be- 
cause they appear along our way of obtaining an effective 
Hamiltonian H (or Ho in the lowest order of the theory) 
in which the statistical degree of freedom <f>i varies in a 
wide range of values (— oo < <pi < oo). This is the only 
consistent interpretation of our consideration performed 
so far. We can assume that up to now we have obtained 
nothing else but an effective free energy (effective lattice 
Hamiltonian) in terms of the new lattice field fa. In a 
clear approximation, this effective model is given by (|6jk 
for the expansion in powers of (pi, see (|T5|) . 

Another important aspect of our consideration is that 
the new Hamiltonian contains more ij— interactions than 
the original model {]]). The mathematical form of (fT5|) 
gives a clear physical interpretation of these interactions: 
the first 4>i4>j— term in the r.h.s. of (JT5J) describes an 
interaction that is quite similar to the original inter-spin 
interaction in ([1]) whereas the second term of the same 
type in (|15p describes an indirect two-site jk— interaction 
that is mediated by the "i— spins" . This means that the 
latter interaction JijJik has twice larger radius of action 
than the original Jy exchange. The four-point interac- 
tion given by the third term in the r.h.s. of (If 5j) can also 
be described in the above style, introduced for the first 
time in this paper. 

Therefore, even at this early stage of consideration 
(H ~ Ho) we see that the effective free energy exhibits ef- 
fects of "statistical extension of the original inter-particle 
correlations (interactions)" . This "principle of growth of 
statistical correlations" is well known in the general phase 



transition theory. Here we show the concrete mechanism 
of the respective phenomenon and a systematic way of de- 
scription of growth of statistical correlations. This point 
will become more clear from the results in the next Sec- 
tion. We shall sec that the investigation along this path 
leads to a quite unexpected and intriguing picture. 
III. BEYOND THE STANDARD THEORY 

A. Perturbation series 

Here we consider the 4>— contributions to the effective 
free energy (Hamiltonian) H(4>) which are generated by 
the term Hf((p). For obvious reasons, terms that are 
(^-independent will be omitted. 

It seems convenient to rearrange our theory by in- 
troducing the auxiliary variables A, = (s^o — (pi), and 
Oi = (si — Sio), where Sio = (si)o- Then H can be writ- 
ten as an infinite perturbation series in powers of the 
(perturbation) Hamiltonian part 

Sf(s, <t>) = JiiWj ~ A «°> ( 26 ) 

ij ij 

The respective series can be presented in the form 

1 oo 

KM = - 2 £ JijAiAj + ]T Hf (</>), (27) 
ij l=l 

where 

^) = ^^(S}( S ,0))oe, (28) 

(...)oc denotes the so-called connected averages [T|; for 
example, the average (Sj)o(Sj) is excluded from the 
connected (Sj)o c - For this cumulant (semi-invariant) ex- 
pansion rules, similar to the Wick theorem in the pertur- 
bation theory of propagator type, are not available and 
one should perform the calculations with some caution. 
The term Wp is equal to zero, and this leads to a re- 
duction of some infinite series in the next orders of the 
theory (I > 1). 

B. Lowest order correction 

Let us consider the first term in the r.h.s. of Eq. (|27|) 
and neglect all others. This yields H = (Ho +Hf) in the 
form 
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1 



n JjjJjkJj 



7~L ~ ^ JjjJik<Pj<Pk 

ijk ijkl 

— — ^ JijJikJilJim4'j ( t ) k4'l < i ) m 
ijklm 



ji9k<Pi 



(29) 



"77" / ^ JijJjkJilJimJii 



»k<Pl<, 



ijkhnn 



Performing this straightforward calculation one readily 
sees a very important property of the present theory, 
namely, that the first term in the r.h.s. of (fTS)) is to- 
tally compensated by a respective counter term coming 
from the new contribution (~ AjAj) to the effective free 
energy. Besides, another new term twice compensates 
the second term in (fT"5|) so that the term of type JJfab 
now appears with a positive sign. The </> 4 — part of the 
effective free energy also undergoes a drastic change due 
to the A A— correction coming from Eq. ([28]) . 

The same result can be obtained in a more general 
and, perhaps, more convenient way, if we add the AA— 
term in (|27p to H.q from © before doing the expansion 
of Landau type. Then, within the same lowest order ap- 
proximation for the series (|27ll we obtain a more general 
result for TL 1 namely 

H = JijtHPa^thipaj)) (30) 

ij 

+ JijfaMPaj) - /T 1 J2 ln [2ch(/9oi)] . 



This form of TL clearly shows the lack of the simple 
Jijfafa term describing the direct two-site exchange. In 
our further considerations we shall be faced only with 
inter-particle interactions (correlations), which are ex- 
tended at distances larger than Ri n t- In order to ob- 
tain (|29|) one must expand the transcendental functions 
in |(3DJ|; [floi) < 1. 

The new forms ([2Uf and (|3"0|) of the effective free energy 
describes only indirect two-site interactions because the 
direct two-site interaction disappeared from our consid- 
eration. Now we are at a stage of description of correla- 
tions which extend up to 2i?^ and larger distances. The 
tendency of the growing length scale of the interactions 
included in our effective free energy will be the main and 



unavoidable feature of our further consideration. 

Another very important feature of our findings is that 
in the simplest variant of the theory when the field is 
uniform (fa = fa) as well as in the GL variant in LWLA 
given by Eq. ([24)) the values of the Landau coefficients Co, 
ro, and uq do not change within the framework of accu- 
racy of the effective 4 -field theory, where terms (Landau 
invariants) of order 0(t 3 ) are ignored; t = (T — T c )/T c . 
This property seems to exist to any order of the expan- 
sion (|27[) in the limit of an infinite-range (z — > oo) initial 
interaction J^. The 1/z— corrections to the parameters 
(co, ro, Mo) of the effective field theory are obtained from 
the (/ > 1)— terms in the series (|27|) . 



C. High-energy corrections of higher order 

We have already mentioned that the (/ = l)-term 
in (|27p is zero. A calculation of the next two terms (/ = 
2,3) in ([2T]) has been carried out in Ref. [2]. Following 

Ref. we can write H^p in the form 



H 



(2) 



j 2 



4 ^ ch^aOch^a,) 
/3 ^ T T AjAfe 



(31) 



ijk 



k 'ch 2 (f3 ai )' 



The result (|3Tj) gives the first (1/z)— corrections to the 
parameters Co, ro, and uq. Adding the result (|3T|) to TL(fa) 
we obtain a form of the effective Hamiltonian TL(fa) which 
is more precise than the preceding ones. Let us write 
down this quite lengthy expression of the effective Hamil- 
tonian: 



7 



02 p3 p3 
H = — JjjJikJjl<f > k<l > l + — JjjJikJil&k&l 77 JijJikJjlJk 



(32) 



ijkl 

4 



ijkl 



ijklm 



^ JijJikJilJim4>3 ( i ) k ( i ) l ( t ) m—f3 A J inJ ijJ ikJ ilJ nnifyj^k^l^m 
ijklm injklm 



^ ^ Jin J ip J 7i j J nkJ nl J prnfyj^k^l^m ^ ^ JinJipJijJikJnlJpm^j^ 



'kVlVrn 



inpjkhn 



inpjklm 



o ^ ] ^in ^itf Jik J il J-i 



im<Pj < Pk(Pl ( Pm 



P 



a ^ ] ^ij ^<fc ^"il ^ 



im9j<Pk9l9m' 



in j kirn 



injklm 



r 



The result (|32| shows that the indirect two-point interac- 
tions of type JJ4>4> available in the effective Hamiltonian 
(f2U]) do not exist in this higher accuracy of the theory. 
The interactions of type JJJ4>(f> and J JJJcjxf) presented 



in the effective Hamiltonian ([32]) extend up to distances 
3Ri n t. The same is valid for the four point interactions 
included in ([52"]) . 

Within LWLA the Eq. ([32]) yields 



J 



0(fci)0(fc 2 )^(fc a )0(-fci - fca - fcs) 



(33) 



r 



In (SI, 



(34) 

where fo = Ub{T — T c ) is given by the "true" (renormal- 
ized) critical temperature 

T c = T c0 (l - ~Y (35) 

In deriving this lattice version of the effective Hamil- 
tonian we have performed the lattice summations in 



J 



in the reciprocal (fe)-space with the help of LWLA: 
J(k) w (J - c k 2 ). 

The present results demonstrate a type of renormal- 
ization of the GL parameters (cq, ro, uq) of the effective 
Hamiltonian due to 1/z-corrections. By a suitable choice 
of units, one of these parameters can be kept invariant, 
for example, equal to unity. Therefore, within a suit- 
able normalization of the theory, the field 4>rx) acquires 
a 1/z-correction as well [2j. 

The term H { p in (27]) has the form 



o/(3) _ P 1 ST T 



th^a^th^aj) 



A z th(/3oj) 



6 2-i ch 2 (/3a i )ch 2 (/3a i )ch 2 (/3a z ) 



Jij J il Jjl 



ijl 



iji 



ck l (Pai)ch''(J3a j 



A, A, 



ijln 



q 7 ^ij u il u ii 



ijln 



A 3 A;A n th(/3a t ) 
ch 2 (/3a,) 



(36) 



r 



(3) 

Let us consider the contribution of the term to the In performing the calculations for both short-range (nn) 
quadratic (<t> 2 -) part of the effective Hamiltonian 
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and long-range (Rmt ^> Go) we have to evaluate again 
several lattice sums. Here we shall mention a particular 
sum, namely, 



j\t Jij J* 1 Jjl ■ 

ijl 



(37) 



which is equal to zero for nn interactions but gives a 
contribution (w J 3 / z) for interaction radius Ri„t > do- 
Thus we introduce the following interpolation formula: 
E = k(z)/z, where < k(z) < 1 is an interpolation 
parameter which is supposed to be a smooth function 
of the coordination number z. The limiting case k = 
corresponds to nn interactions and the limiting case k — 
1 corresponds to interactions of larger size. We suppose 
that the shape of the function k(z) depends on details of 
the function J(R). 

Bearing in mind these notes we have calculated the 
quadratic (jfiiifij—) contribution to the effective Hamil- 
tonian Tl(4>) which comes from the Hamiltonian part 
Hf((j)) given by (f36|) . Here we present the results for 
the parameters T c q, c, and r which define the quadratic 
part H 2 of H: 



n 2 (<f>) 



(38) 



k 



with the (bare) correlation function 

Gu\k) =X>i(z) [T-T c +V 2 (z)T cP 2 k 2 ] . (39) 

Here 

1 + k 1 + 3k\ „ 

t,.= \i-— z ^r) T -^ ( 4 °) 



V x (z) = 1 + 



3 + 4k 22 + 60k + 30k 2 



and 



V 2 {z) = 1 



3z 2 



2 + 3k 14 + 30k + 9k 2 



3z 2 



(41) 



(42) 



The result ([55)1 for T c (z) has been published for the first 
time in case of nn interactions (k = 0); note, that 
there are errors in Ref. [3j for the functions Z?i(z) and 
V 2 (z). 

The functions T>\(z) and T> 2 (z) renormalize the field 
4>(k) and the vertices Co, and vq. For a total renormal- 
ization of the parameters of the theory up to the second 
order in the 1/z— expansion we need to know the (1/z) 2 - 
correction to the vertex uq- We suppose that the calcu- 
lation of this correction can be accomplished in the way 
described above; this "z-renormalization" has been dis- 
cussed to first order in (1/z) in Ref. 0. Here we wish to 
stress that within our extension of the theory the mag- 
netic susceptibility Go(0) is T>\ times smaller than the 
known MF susceptibility corresponding to 2?i(oo) = 1. 



The numerical coefficients in (j40|) - (|42|) indicate that 
real numbers z of nn like n = 2, 4, 6 for simple lattices 
of spatial dimensionalities D — 1,2, 3, respectively, give 
a good expansion parameter 1/z. The 1/z-orrections are 
more substantial for the case of short-range interactions 
(Hint ~ a o): an d one may suppose that for such inter- 
actions the (l/z)-series (|4U| - (l4"2"j) are asymptotic; for 
the case of T c , see a discussion of this topic in Ref. 0]. 
But even in case of asymptotic type of these series they 
may give more reliable results than the "bare" values 
(co, 7*o, uq) of the Landau parameters; see arguments pre- 
sented in Ref. M- 



D. Critical temperature 

The critical temperature T c given by (|40[) can be 
compared with exact and reliable numerical (MC) re- 
sults. Let us consider nn interactions (k = 0). For 
one-dimensional (ID) IM, we know that T c = 0, MFA 
predicts T c = 2Jo/fcs (in this case, z — 2D = 2), 
and ([4*0]) yields T c — 5Jo/6fcs. This is a quite good 
result for ID systems with very strong fluctuation ef- 
fects. In 2D systems the fluctuations are not so strong 
and we find that (|4"0"]) reproduces the exact Onsager re- 
sult (T c — 2.27Jo/fcs) with an error of 22%, i.e., we 
have T c = 35Jo/12fcs. For 3D systems our result is 
T c = 89Jo/18fcs, whereas the best series analysis and 
MC results yield a difference of 9%: T c = 4.5Jo/fcs (see, 
e.g., Refs.[l8l [l9l]). Our results seem quite reliable. Let 
us emphasize, that the 1/z-series, as almost all most rele- 
vant series, known in theoretical physics, is an asymptotic 
series. Therefore, one may expect, that the results for T c 
will be worsened after some order in 1/z, for example, 
fourth order for T C (3D), and third, or, even second order 
for T C (2D). 



E. Ground state 

The l/z 2 -correction to the vertex uq has not been 
yet calculated although this calculation does not present 
difficulties. In this situation we shall give notion for 
the ground state energy by using the first order cor- 
rections to ro and uq. The equilibrium free energy 
per site / = (F/N) is given by F = H and j33]) 
for the (k — 0)-Fourier amplitude 0(0), which mini- 
mizes /. Denote for convenience <f>(0) = \/~Nip ^ 
for the low-temperature ordered phase. From (33) we 
obtain / = — (r 2 /16w) < whereas the "unrenormal- 
ized free energy is fa = (Ho/N) = —(r 2 /16uo) < 
Thus, using J33) we have f(T) = (1 + 2/z)f (T), which 
means that the new effective theory has a lower energy 
of the ordered phase. This is true also in the case of 
T = 0, where tq = —UbTcq = — J. For the zero tem- 
perature (ground) state we have /o(0) = — 3J/4 and 
/(0) = -3(1 + 2/z) J/4. This result is also along the 
right direction because the MF theories (7io) give unre- 
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liable high values of the ground state energy. The order 
parameter (fi 2 (T) = (— r/4u) is (1 — z) times smaller than 
the respective quantity <fo(T) = (— r /4u ) in the usual 
theory based on Ho- 



F. Effective interactions and growing of fluctuation 
correlations 

We have explicitly shown the phenomenon of the grow- 
ing length size of the interpaticle correlations in a classic 
system of interacting particles. To see this we have al- 
ready introduced a new interpretation of the terms in 
the effective Hamiltonian (see Sec. II. D). Let us consider 
the terms present in 7i as terms describing certain inter- 
site interactions. While the initial interaction Jy in IM 
ensures only two-site correlations (interactions), the ef- 
fective Hamiltonians P^)l . (f2"9")) , and (|3"!2)) contain multi- 
site effective interactions. In contrast to the usual theory 
(|15[) . where only extremely short-range effective corre- 
lations are contained, the more precise effective Hamil- 
tonians contain long-range two-site (4>i4>j) and four-site 
(4>i(j)j(f>k4 l i) correlations, and all of these correlations are 
indirect, i.e. the correlation, for example, between two 
sites (ij) is mediated by one or more other sites (fc, ...). 
A direct ( )-interaction is presented by the first term in 
the r.h.s. of (fT5f but also the system exhibits two indi- 
rect correlations of type <j)i<j)j and <pi4>j4>k4'i given by the 
last two terms in the r.h.s. of (|15p. In the more precise 
variants of the theory, where a larger portion of the initial 
partition function has been calculated, the direct inter- 
site interaction vanishes, and the particles are correlated 
only by indirect effective interactions. The length scale 
of these correlations grows in a monotonous way with the 
increase of the accuracy of the calculation, i.e. with the 
increase of the number I of the terms in the series (|27|l . 
If we take the two-site correlations in the nn IM as an 
example, the maximal length of extension of these corre- 
lations in |TJ is oto, in ^ - 2a , in (JUJ) - 3a , in ^ - 
4ao, i.e. (p— l)ao, where p is the maximal number of the 
summation indices in the terms of type (fri'fij in a given 
effective Hamiltonian. Surely, the number p tends to N. 
This means that the most accurate effective field theory 
of many-body systems will correspond to (almost-) infi- 
nite range of correlations. 

The origin of these correlations is purely statistical. 
This effect is known and has both general formulation 
and application in many-body physics. Here we have 
established and described in details the concrete mech- 
anism of this effect and, moreover, we have performed 
a demonstration of the remarkable picture of successive 
growth of the correlation length scale. 



G. Final remarks 

Obviously the (l/z)-corrections are not the main point 
of discussion at the end of this paper of restricted 



length. Let us mention that the growth of the corre- 
lations discussed in Sec. III.F is not related with the 
l/z-corrections. It exists for any z, even in the "MF lim- 
iting case" of z — > oo, when the GL parameters keep their 
"initial" values (co, ro, uo)- This effect follows from the 
fact that the terms in the initial Hamiltonian are com- 
pensated by the first "fluctuation" correction; see the first 
term in the r.h.s. of (|27p. At the next level of accuracy 
of the calculation, terms coming from the (I = 2)-term 
in l|27p compensate the available terms and this process 
continues up to the incorporation of all particles in the 
correlation phenomenon; remember that the term corre- 
sponding to I = 1 is equal to zero. 

Here we emphasize that the terms in H.q - actually 
one of the most often used Hamiltonian, does not exist 
at all. They vanish just after the inclusion of the first cor- 
rection to the usual theory; see AjA^-correction in l|27[). 
In place of these terms other terms with more complex 
structure come from the perturbation series l|27p. The 
outlined picture clearly indicates, that the terms which 
finally remain in the <^ 4 -theory, are terms of type 

- {I3 M - 1 J M - )3 M J M+1 ) <j?, M~N, (43) 

where obvious notations have been introduced; for exam- 
ple, (3 1 J 2 denotes the first term in the r.h.s. of ([2"9"p . The 
</> terms behave differently, because a lowest order term 
in J, namely, a term of type (3 3 J^cf) 4 appear at any step 
of development of the series (j27|) . 

At any stage of this surprising picture of the infinite 
series of successive modifications of the Hamiltonian both 
(f) 2 - and c/) 4 -terms keep their numerical coefficients equal 
to that in the usual GL Hamiltonian TCq. This is true 
within the whole scope of validity of the expansion in 
powers of <p. 

An important note, which should be emphasized is the 
following. While the sum ([2]) is invariant with respect 
to the site i in regular lattices, the sum ([7|) depends on 
the site i. The reason is that the field configuration {4>i\ 
which takes part in Eq.([7]) is not the equilibrium field. 
For the equilibrium field cf>i the sum ([7]) will not depen- 
dent on the site i. This is consistent with the general 
notion that the equilibrium order in the volume of a ho- 
mogeneous system in lack of effects of external fields, 
should be uniform. 

Our consideration justifies the GL fluctuation Hamil- 
tonian. However, we have presented a new and quite 
surprising picture of the inter-particle correlations, 
which reveals new remarkable properties of the GL 
theory. Apart from the 1/z-corrections to the GL 
parameters, the structure of this theory is absolutely 
comprehensive as a tool for investigation of large-scale 
correlation phenomena in many-body systems. We are 
certain that our findings have an application beyond the 
field of phase transitions. 
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